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Two-electron correlated spectra of non-sequential double ionization below laser-intensity thresh- 
old are known to exhibit back-to-back scattering of the electrons, viz., the anticorrelation of the 
electrons. Currently, the widely accepted interpretation of the anticorrelation is recollision-induced 
excitation of the ion plus subsequent field ionization of the second electron. We argue that another 
mechanism, namely simultaneous electron emission, when the time of return of the rescattered elec- 
tron is equal to the time of liberation of the bounded electron (the ion has no time for excitation), can 
also explain the anticorrelation of the electrons in the deep below laser-intensity threshold regime. 
Our conclusion is based on the results of the numerical solution of the time-dependent Schrodinger 
equation for a model system of two one-dimensional electrons as well as an adiabatic analytic model 
that allows for a closed-form solution. 
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I. INTRODUCTION 

After one-electron ionization of an atom or a molecule, 
the free electron can recollide with its parent ion [1, 2], 
Inelastic scattering on the ion may result in collisional 
excitation or ionization. Liberation of the second elec- 
tron during recollision is known as non-sequential double 
ionization (NSDI). 

Originally, NSDI has been discovered experimentally 
for alkaline-earth atoms [3] (regarding further experimen- 
tal investigations of NSDI for alkaline- earth atoms, see, 
e.g., Refs. [4-9] and references therein). Later, NSDI has 
been found in nobel gas atoms [10, 11]. 

According to classical calculations, the maximum ki- 
netic energy of the recolliding electron is sw 3.2t/ p [2], 
where U p = (F/2lu) 2 , F is the laser field strength, and 
lu is the laser frequency [unless stated otherwise, the 
Hartree atomic units (a.u., h = m e — |e| = 1) are used 
throughout the paper]. Hence, one can divide NSDI into 
two types: when the energy of the recolliding electron is 
enough to collisionally ionize the parent ion - the above 
laser-intensity threshold regime, and when the kinetic en- 
ergy is insufficient to directly ionize the ion - the be- 
low laser-intensity threshold (BIT) regime. The former 
one regime is thoroughly studied experimentally as well 
as theoretically (see, e.g., Refs. [12-30] and references 
therein). 

Despite tremendous challenges, NSDI BIT has been 
experimentally observed [31-36]. The main feature of 
this regime, revealed by measured correlated two-electron 
spectra, is that the two electrons prefer to drift out in 
opposite directions, manifesting the so-called anticorre- 



lation of the electrons. The interpretation of this ob- 
servation is the aim of recent theoretical studies [37-43] . 
Briefly summarizing these models, it is widely accepted 
that recollision-induced excitation of the ion plus sub- 
sequent field ionization of the second electron (RESI) 
is a main mechanism of NSDI BIT. Nonetheless, alter- 
native, non-RESI, mechanisms are also considered (see, 
e.g., Refs. [30, 35, 42, 44] and references therein). 



'Electronic address: dbondarSsciborg .uwaterloo . ca 
^Electronic address: gennady . yudinOnrc . ca 



However, is excitation of the parent ion indeed nec- 
essary to explain the electron anticorrelation in NSDI 
BIT? We show that this is not always the case. Our con- 
clusions are based, first, on model ab-initio calculations 
for a system of two one-dimensional electrons show that 
the anticorrelation of the electrons exists even if the ion 
has only a single bound state (Sec. II). Second, within 
the adiabatic approximation (see, e.g., Refs. [45-53] and 
references therein), we present a simple analytical three- 
dimensional model based on the assumption that both 
the electrons are ejected simultaneously (the time of re- 
turn of the first electron coincides with the time of libera- 
tion of the second electron, i.e., the ion has no time to be 
excited) in Sec. III. An advantage of this model is that it 
allows for a simple analytical solution in closed form. In a 
certain range of parameters, the correlated two-electron 
spectrum obtained within this model exhibits the anti- 
correlation of the electrons. This mechanism of simulta- 
neous electron emission can produce the anticorrelation 
of the electrons in the deep BIT regime. 
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II. AB INITO EVIDENCE OF THE NEW 
MECHANISM 

Consider the model Hamiltonian for a system of two 
one-dimensional electrons 



H(t) 



(Pi 
-V(x 2 



it) 



/2 + V(xi) + V(x 2 ) 
Xx) + (xi + x 2 )F L (t), 



(1) 



where V(x) — — 4exp(— 3a; 2 ) is the prototype for the 
potential of the electron-core attraction, — V(x 2 — x%) 
is the prototype of the electron-electron repulsion, and 
F L (t) = Ff{t) sin(wt) (F = 0.05 and w = 0.6) - the laser 
pulse, where f(t) is a trapezoid with one-cycle turn-on, 
six-cycle full strength, and one-cycle turn-off. The po- 
tential V(x) is chosen such that the one-particle Hamil- 
tonian, Hi on = —d 2 /(2dx 2 ) + V(x), supports only one 
bound state with the ionization potential I p = 2.11. 

To make analysis more transparent, we substitute the 
original problem (NSDI) by the corresponding problem 
of laser-assisted scattering, i.e., we simply discard the 
first step of NSDI - liberation of the first electron. In 
other words, instead of assuming that the two-electron 
system initially is in the ground state, we assume that 
the first electron is an incident wave-packet and the sec- 
ond electron is in the single bound state of the ion. This 
modification is very useful, as it allows us to eliminate all 
other possible interactions and processes except the three 
major components - the electron-electron repulsion, the 
electron-core attraction, and the electron-laser field in- 
teraction. 

We solve numerically, by means of the split-operator 
method, the time-dependent Schrddinger equation, 

id^(xx,x 2 ;t)/dt = H{t)^(x u x 2 ;t). (2) 

The initial unsymmetrized wave function takes the form 

(Si - [if 



^(xi, x 2 ; 0) = exp 



(2a x 



i\/2E v , 



Xi 



1>(X2), (3) 



where a x = 2, fi = 7a x , H ion tp(x) = -I p ij)(x), and 
Ei n = 0, 0.3 is the mean value of the kinetic energy of 
the incident electron. The parameters are selected such 
that, I p -(F/u) 2 /2 « 4oj and/ p -(V(16+F/a;) 2 /2 ~ 3w, 
i.e., the bound electron needs to absorb at least four 
(three) photons to be liberated in the case of E in = 
{Ei n = 0.3). The wave function then is symmetrized, 
i.e., we assume that the spins of the electrons are an- 
tiparallel. 

The wave function of the model system is pictured in 
Fig. 1 for two values of Ei„. The wave function in the 
absence of the laser produces only maxima on the axes 
xi.2 corresponding to one electron bound and the other 
one free; hence, this part of the wave function, which also 
shows up when the laser field is turned on, should be ig- 
nored because it does not correspond to NSDI. When the 
mean kinetic energy of the incident particle is zero, Fig. 
1(a), we observe that the two electrons "prefer" to be 
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FIG. 1: (Color online) Plots of the wave functions of the 
model system, log \^(xi, X2', 73) | 2 , for different values of the 
kinetic energy of the incident electron [Eqs. (1) and (2)]. (a) 
E m = 0; (b) E in = 0.3. 



anticorrelated rather than correlated. However, once the 
mean kinetic energy is slightly increased such that NSDI 
still progresses in the BIT regime [see Fig. 1(b)], the elec- 
tron correlation begins to dominate the anticorrelation. 
Note that the wave function in the third quadrant of Fig. 
1(b) (i.e., the portion of the wave function that describes 
the correlation of the electrons) has a characteristic V- 
like (aka fingerlike) shape, which is a dominant shape of 
the two-electron correlated spectra of NSDI in the above 
threshold regime. 
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III. THE SIMULTANEOUS ELECTRON 
EMISSION WITHIN THE ADIABATIC 
APPROXIMATION 

We calculate the correlated spectrum within the 
adapted Savichev adiabatic approximation [48, 51] 



/<Po 
1 



r oc cxp Im / [Ef((p) — Ei(ip)]dip 



(4) 



which is a generalization of the Dykhne adiabatic method 
[45] to the case when a quantum system has a continuum 
spectrum. Here, p = cot is the phase of a laser field, ifo 
is the complex solution of the equation Ei(ip ) = Ef(<po) 
with the smallest positive imaginary part, Eij-(p) being 
the adiabatic terms of the system "before" and "after" 
a non-adiabatic transition. If the quasi classical approx- 
imation to T is sufficient, then instead of using the adi- 
abatic terms, one can employ the total energies of the 
corresponding classical systems [51]. 

As previously, we replace the problem of NSDI BIT 
by the problem of laser-assisted scattering of an electron 
by an ion. We study the simultaneous electron emission 
(SEE) process, when the moment of collision of the inci- 
dent electron coincides with the moment of ionization of 
the ion. 

To include electron-electron repulsion, we apply the 
standard exponential perturbation theory. Namely, first 
we find electron energies and trajectories without the 
electron-electron repulsion. Then, we correct electron 
action and energies by adding the effect of the electron- 
electron repulsion, which calculated along the zero-order 
trajectories. 

Thus, in zero order, we define Eij(tp) without the 
electron-electron repulsion. Since before ionization, one 
electron is free and the other is bound, the total classical 
energy of the system before collision is 



EM = [p + A(^)] 2 /2-/ P 



(5) 



J 



where p is the (three-dimensional) canonical momentum 
(i.e., the kinetic momenta at p = ±oo) of the incident 
electron, I p is the ionization potential of the ion, and 
A(ip) = — (F/w) simp is the vector potential of a linearly 
polarized laser field. After collision both the electrons 
are free, and the classical energy of the system reads 



E f (p) = [k 1+ A^)] 2 /2+[k 2 + A(^)] 2 /2, (6) 



where ki j2 are canonical momenta of the first and second 
electrons. 

Such two electron process is formally equivalent to 
single-electron strong field ionization of a quasiatom 
(within a pre-exponential accuracy). This statement is 
manifested by the following equality 



I r<po 
u J 



[Ef(<f) ~ Ei(ip)] d<p 



^°|I[jr + ^)] 2 + ^J#, (7) 

where the right hand side of Eq. (7) is the action of 
single-electron ionization within the strong field approx- 
imation, — k\\ 1 + k\\ 2 —p\\ is the effective longitudinal 
momentum, and J v = I p + (k 2 + k 2 . p 2 Jf 2 ) /2 is 
the effective ionization potential of the quasiatom. (Note 
that such a reduction is in fact more general and can be 
applied to an arbitrary number of electrons, for further 
discussions see Ch. 7 of Ref. [54].) Therefore, derivation 
of the correlated spectra for SEE is reduced to calcula- 
tion of the momentum distribution of photoelectrons af- 
ter single-electron ionization without any restrictions on 
Jif . The most suitable solution of the last problem for 
our current discussion has been found in Ref. [55] . Note 
that the equation obtained in Ref. [55] is valid only for 
positive J v ; hence, the case of J p < being of interest 
for SEE, must also be considered. 



Substituting Eqs. (5) and (6) into Eq. (4) and taking into account the previous comments, we obtain 

r(ki, k 2 ) ex exp (-2 \S P \ G/u + V ee ) , 

J arcsinx_ +«arccoshx + if ,J? p < 
^° 1 arcsiny_ + iarccosh y + if ,_^ p > 0, 

— 1 ) arccoshx + 




' x \-l{ 


' 1t)X- 


f x + 


l-2x_ 






ly\ i 




+ y+ 


l-2y 2 _ 




2s 2 



2p- + lj arccoshy + 

*± = |g(7?-l) + l| /2 ± |g (77 - 1) - 1| /2, 
V± = V(OT + l) 2 +.9 2 /2 ± V(OT-l) 2 +5 2 /2, 



if J v < 0, 
) if^ P >0, 



(8) 
(9) 

(10) 



g = Lo^2\Jp\/F 7 V = jr/v/2pg. 

In Eq. (8), V ee denotes a crucial correction due to the electron-electron repulsion, 

^° dp _ -2 T _ f va dp 

|ri(^)-r 2 (^)| 



V e , 



= -iim r 

u Jr ( 



>o k x 



-Im 



JRc V „ y/(ip ~ p ) 2 +£ 2 



-it/ |ki 



(11) 



k ¥ (a.u.) 



fcp (a.u.) 
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FIG. 2: (Color online) Correlated two-electron spectra of Ar 
(linear scale) given by Eq. (8) at 800 nm and 1 x 10 13 W/cm 2 
(p_i_ = k±i = k± 2 — 0) for different momenta of the incident 
electron: (a) py = 0.1, (b) py = 0.2, (c) py = 0.25, (d) py = 
0.4. 

The contribution of the Coulomb potential has been 
calculated along the field- free trajectories [39, 56] 

i r 

ri,2(y) = - / [ki, 2 + M<t>)] d<t>- 

The condition of validity of Eq. (8) is T <C 1, which 
is merely the condition of applicability of the adiabatic 
approximation. 

The model of NSDI BIT presented here [Eq. (8)] is 
similar to the one reported in Ref. [39]. However, there 
is an important difference. In the correlated spectra (8), 
the momentum of the incident electron p is a free pa- 
rameter, whereas the correlated spectrum given by Eq. 
(19) of Ref. [39] does not have this freedom - the canon- 
ical momentum of the recolliding electron is fixed as a 
function of the phase of recollision (this functional de- 
pendence is obtained from the saddle-point analysis of 
the S-matrix amplitude for NSDI). Therefore, Eq. (8) 
allows us to establish links between values of the canoni- 
cal momentum of the rescattered (incident) electron and 
specific portions of the correlated spectra. The depen- 
dance of the correlated spectra (2) on the momentum of 
the incident electron is presented in Fig. 2. 

The peaks of the momentum distributions depicted in 
Figs. 2(c) and 2(d) are located approximately at 10^/Up. 
This value is much larger than 2 \fUp that is the classical 
drift momentum of a free electron in a laser field. Such 
a discrepancy is due to the fact that the momentum of 
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FIG. 3: (Color online) Correlated two-electron spectra of Ar 
(linear scale) given by Eq. (8) at 800 nm (p = and k± 1 — 
k±2 =0) for different intensities of the laser field: (a) 3 x 10 13 
W/cm 2 , (b) 2.7 x 10 13 W/cm 2 , (c) 2.5 x 10 13 W/cm 2 , (d) 
2.3 x 10 13 W/cm 2 , (e) 2 x 10 13 W/cm 2 , and (f) 1 x 10 13 
W/cm 2 . 



the incident electron p can take an arbitrary value in our 
model, and its value is not determined by the saddle- 
point integration of the corresponding S-matrix ampli- 
tude, which is in turn equivalent to the employment of 
the quasi-classical description of the rescattered electron. 
(As we have mentioned, the incident electron momentum 
was fixed in Ref. [39], and hence, such high electron mo- 
menta were not observed.) 

Comparing Figs. 1 and 2, one concludes that the ana- 
lytical model [Eq. (8)] qualitatively agrees with the nu- 
merical solution of the Schrodinger equation (2) for the 
model system of two one-dimensional electrons presented 
in Sec. II. Indeed, when the momentum of the incident 
electron is zero - the anticorrelation of the electrons dom- 
inates the correlation [Figs. 2(a)]. The situation is op- 
posite when the momentum of the incident electron is 
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slightly increased [Figs. 2(c) and 2(d)]. 

An interesting question is the dependence of the cor- 
related spectrum (8) on the intensity of the laser field, 
which is pictured in Fig. 3. The parameters used to 
plot Fig. 3(a) coincide with the parameters employed 
in the recent experiment [36]. Yet, the experimentally 
measured correlated spectrum exhibits peaks in the sec- 
ond and fourth quadrants, which contradicts Fig. 3(a). 
The reason of such a disagreement is that the anticor- 
relation in these experimental data is due to the RESI 
mechanism. As the intensity lowers, the peaks in the 
correlated spectra of SEE shift to the second and fourth 
quadrants [see Fig. 3(f)], i.e., the SEE process leeds 
to the anticorrelation of the electrons in the deep BIT 
regime. This can be explain intuitively in the following 
way: The lower the intensity, the smaller the canonical 
momentum of the returning (incident) electron. Since 
the canonical momentum of the system is approximately 
conserved, we obtain w ki + k 2 ; hence, ki w — k 2 . 
Having absorbed a necessary number of photons, both 
the electrons emerge in the continuum where they ex- 
perience strong electron-electron repulsion that pushes 
them apart. Therefore, during SEE, the electrons gain 
momenta because of the electron-electron repulsion. In- 
deed, if the electron-electron interaction is "turned off" 
in Eq. (8) by setting V ee = 0, then we obtain a single 
peak centred at the origin instead of the two peaks visible 
in Fig. 3(f). Recall that contrary to SEE, the electrons 
gain kinetic energy by oscillating in the laser held in the 
course of RESI. 



IV. CONCLUSIONS AND DISCUSSIONS 

We have demonstrated that the mechanism of simul- 
taneous electron emission, when the time of return of the 
rescattered electron is equal to the time of liberation of 
the bound electron, can be responsible for the anticorrela- 
tion of the electrons during NSDI in the deep BIT regime 
[see Fig. 3(f)]. The SEE process significantly differs from 
RESI because it does not require the excitation of the ion 
to explain the anticorrelation of the electrons observed in 
the two-electron correlated spectra. According to SEE, 
NSDI BIT can be realized by a rescattered electron with 



a very low (canonical) momentum. Such an electron lo- 
cated near the parent ion creates strong Coulomb field 
that acts on the bound electron. This Coulomb repul- 
sion eventually ionizes the parent ion and leads to the 
anticorrelation of both the electrons. The SEE should 
be manifested at small laser intensities because the lower 
the strength of the laser field, the lower the canonical 
momentum of the rescattered electron. 

Another important difference between SEE and RESI 
mechanisms is that the former results strictly in the an- 
ticorrelation of the electrons, viz., the correlated two- 
electron spectra should be concentrated mainly in the 
second and fourth quadrants. Whereas, the RESI mech- 
anism leads to electron momentum distributions more or 
less equally distributed in all four quadrants of the mo- 
mentum plane [40, 41]. 

The main difference between SEE and other proposed 
non-RESI mechanisms (see, e.g., Refs. [35, 42, 44]) is 
that it is induced by rescattered electrons with nearly 
zero momentum, and hence it requires no time delay be- 
tween the return of the first electron and liberation of the 
second. 

SEE and RESI are by no means mutually exclusive 
processes of NSDI; they both contribute. It is important 
to study quantitatively the relative contribution of SEE 
and RESI to NSDI in the deep BIT regime. Nevertheless, 
this cannot be done by neither the reduced dimensional 
numerical model presented in Sec. II, nor by the analyt- 
ical adiabatic model developed in Sec. III. 

Our analytical model can be further developed in the 
spirit of the correlated Keldysh-Faisal-Reiss theory [57- 
59]. However, if one wants to study NSDI in very weak 
laser fields, the model may need a substantial revision 
because in such a regime, the Coulomb interaction dom- 
inants over the interaction with the laser, as a result, 
the classical trajectories of the electrons after ionization 
become chaotic (see, e.g., Refs. [60, 61]). 
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